APS/123-QED 



Width of Percolation Transition in Complex Networks 

Tomer Kalisky^'Q and Reuven Cohen^ 

^Minerva Center and Department of Physics, 
Bar-Ilan University, 52900 Ramat-Gan, Israel 
^Department of Computer Science and Applied Mathematics, 
Weizmann Institute of Science, Rehovot, Israel 
(Dated: February 2, 2008) 

Abstract 

It is known that the critical probabihty for the percolation transition is not a sharp threshold, 
actually it is a region of non-zero width Apc for systems of finite size. Here we present evidence 
that for complex networks Apc ~ where / ~ N'^°p* is the average length of the percolation 
cluster, and N is the number of nodes in the network. For Erdos-Renyi (ER) graphs Uopt = 1/3, 
while for scale-free (SF) networks with a degree distribution P{k) ~ fc^^ and 3 < A < 4, Uopt = 
(A — 3)/(A — 1). We show analytically and numerically that the survivability S{p,l), which is the 
probability of a cluster to survive I chemical shells at probability p, behaves near criticality as 
S{p,l) = S{pc,l) • exp[{p — Pc)l/Pc\- Thus for probabilities inside the region \p — Pc\ < Pc/l the 
behavior of the system is indistinguishable from that of the critical point. 
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I. INTRODUCTION: 



Recently the subject of networks has received much attention. It was reahzed that many 
systems in the real world, such as the Internet, can be successfully modeled as networks. 
Other examples include social networks such as the web of social contacts, and biological 
networks such as the protein interaction network and metabolic networks . The prob- 

lem of percolation on networks has also been studied extensively (e.g. IJ])- Using percolation 
theory we can describe the resilience of the network to breakdown of sites or links ^, 0] , 
epidemic spreading 9], and the properties of the optimal path in a highly network with 
highly fluctuating weights on the links ^. 

A typical percolation system consists of a d-dimensional grid of length L, in which the 
nodes or links are removed with some probability 1 — p, or are considered "conducting" 
with probability p (e.g. P, Below some critical probability Pc the system becomes 

disconnected into small clusters, i.e., it becomes impossible to cross from one side of the grid 
to the other by following the conducting links. Percolation is considered a geometrical phase 
transition exhibiting universality, critical exponents, upper critical dimension at = 6 etc. 
It was noted by Conigilio that for systems of finite size L the transition from connected 
to disconnected state has a width Apc ~ jj/i^y where u is a. critical exponent related to the 
correlation length. 

Pe..co.atio„ o. .etwo... was .tud.ed also .o-n a mathematical po„. of view fl y Q. 

It was found that in Erdos-Renyi (ER) graphs with an average degree (k) the percolation 
threshold is: Pc = Below Pc the graph is composed of small clusters (most of them 
trees). As p approaches pc trees of increasing order appear. At p = pc a. giant component 
emerges and loops of all orders abruptly appear. Nevertheless, for graphs of finite size N it 
was found that the percolation threshold has a finite width Ap^ ~ ^^173 JilJ, meaning that 
all attributes of criticality are present in the system in the range [pc — ^PcPc + ApJ- For 
example: The number of loops is negligible below Pc + Ap^ [s^. 

In this paper we study the Survivability of the network near the critical threshold. The 
survivability S{p,l) is defined to be the probability of a connected cluster to "survive" up to I 
chemical shells in a system with conductance probability p IJ] (i.e the probability that there 
exists at least one node at chemical distance I from a randomly chosen node on the same 
cluster). At the critical point Pc, the survivability decays as a power-law: S{pc,l) ~ 
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where x is a universal exponent. Below pc the survivability decays exponentially to zero, 
while above pc it decays (exponentially) to a constant. Here we will derive analytically and 
numerically the functional form of the survivability above and below the critical point. We 
will show that near the critical point S{p,l) = S{pc,l) ■ exp[{p — Pc)^/Pc\- Thus, given a 
system with a maximal chemical length /, for probabilities inside the range \p — Pc\ < ^ 
the behavior of the system is indistinguishable from that of the critical point. Hence we get 
Ap, ~ s^. 

The maximal chemical length / at criticality is actually the length of the percolation 
cluster, which was found to be: / ~ ]\f'^°pt where is the number of nodes in the network. 
For Erdos-Renyi (ER) graphs Uopt = 1/3, while for scale- free (SF) networks with a degree 
distribution P{k) ~ k'^ and 3 < A < 4, Uopt = (A - 3)/(A - 1) 3. 



II. ERDOS-RENYI GRAPHS: 



Consider an ER graph with a mean degree {k) , and each link having a probability p to 
conduct. We define Ni{x) = uq + riix + n2X^ + n^x^ + ... to be the generating function of 
the number of sites that exists on layer (i.e. chemical shell) / starting from a random node 
on the graph (for a conduction probability p) . 

The generating function for the degree distribution of a randomly chosen node in the 
network is Gq{x) = ^ P{k)-x^ and the generating function for the number of links emerging 
from a node reached by following a randomly chosen link is Gi{x) = T.^kP{k) ■ x^~^ [l5[- 
Taking into account the probability p for conduction, we have: 

Giix) = l-p + py2lTTkPik)-x''-\ (1) 

We can now write the following recursive relation \\^ : 

Ni+,{x) = G,iNi{x)), (2) 

which means that the probability nf'^^'^ for having i nodes at layer / + 1 is composed of 
the probability of reaching a vertex by following a link, and reaching i nodes at layer / by 
following all branches emerging from that vertex - see sketch in Fig. Q 

It can be seen that Ni{0) = uq is the probability that there are nodes at layer /, i.e., 
the probability to die before layer /. Thus e/ = 1 — A^/(0) is the probability to survive up to 
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layer /. From (j2I) we have: 



iV;+i(0) = Gi(iVz(0)) 



(3) 



=Gi(1-q) (4) 

For ER graphs ^0(0;) = Gi{x) = e<^>(^-^) (for p = 1), Thus: 

1 - ei+i = l-p + pe<^>(^-^'~^) = l-p + pe-^''^" (6) 



Q+i =p-pe 



(7) 



Setting 6 = p — Pc, where Pc = and expanding by series we get: 



ei+i = pe + 5 - (Pc + 5) ( 1 - {k)ei + -{k)h} 
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(8) 



Leaving only expressions up to second order in 5 and (we assume that p < Pc and thus 
<^ 1 for large /) we get: 

ei+i^e,-l{k)ef + 6{k)ei (9) 



dei 1 , , V 9 5 

At criticality, 5 = and the solution to this equation is: e; ^ 
suggests the following solution near criticality: ei ~ ■ exp 
In terms of survivability this can be written as: 

S{p,l) = S{pc,l) - exp ( — {p-Pc)l 



(10) 

The additional term 



^6l) 2M 

Pc J ' ' 



(11) 



In order to check this result we numerically solved the survivability S{p, I) near pc according 
to the exact enumeration method presented in ^ 2^. Fig. |21 shows the survivability S'(p, /) 
for different values of p. For p = Pc the survivability decays as a power law, while above 
and below there is an exponential decay, either to zero (for p < pc) or to a constant (for 
p > Pc)- Fig. 01 shows that all curves of the survivability S{p, I) from Fig. El can be rescaled 
such that they all collapse. Moreover, scaled survivabilities from all different graphs with 
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different values of {k) (i.e., different values of pc) also collapse on the same curve. However, 
equation (fTTj) is true only below the percolation threshold where there is no giant component. 
Above the percolation threshold there is an exponential decay to a non-zero constant, and 
the generalized expression is: 

S{p,l) = S{pc,l) -exp ( - — \p - Pc\n +pPoo, (12) 



\ Pc / 

Where pPoo is the probability for a randomly chosen site to be inside the percolation clus- 
ter |23|. Indeed, setting e/+i = ei in equation (|Zj) the resulting "steady state" solution is 



pPoc mm- 



III. SCALE-FREE GRAPHS: 

Scale-free graphs can^be taken to have a degree distribution of the form P{k) = ck~'^ 
where c ~ (A — l)m^^^ [6]. In order to solve equation (0) we have to evaluate: 

G,{l-e,) = ^^Y.^Pmi-e,f-' (13) 

Expanding by powers of e;, and inserting P{k) = ck~^ with 3 < A < 4, we get Q]: 

J2 kP{k){l - ei)'-' ^ {k) - {kik - l))ei + |r(4 - \)ef-' (14) 

Thus equation ^ becomes: 

1 - e,+i ^ 1 - P + ^ {{k) - {k{k - 1))Q + |r(4 - A)e^^) (15) 

Taking p = Pc + S: 

1 - Q+i ^l-{p. + S) + ^ ((A:) - {k{k - 1))Q + ^r(4 - A)ef-2) . (16) 
Substituting pc = (k{k-i)) B S^^- 

1 - « 1 + P=4r(4 - A) . - ^ • + 4r(4 - A) . (17) 

Setting A = pc^T{4: - A) we get: 



Q+i -ei^-A- + --6ei--- 6ef-^ (18) 

Pc Pc 
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Pc 



(19) 



For large /, e; <^ 1, and taking into account that A — 2 > 1 we have ^ ^ e;. Therefore: 



dl 



Pc 



6e,. 



(20) 



. The additional term suggests the 
2511. A similar form can be found for 



For 6 = the solution is ei ~ with x = |l6 
following solution near criticality: e/ ~ ■ exp (^f-^ 
A > 4 j26|. The scaling form for SF networks is also confirmed by numerical simulations as 
shown in Figures |3 and El 



IV. SUMMARY AND CONCLUSIONS 



The scaling form of the survivability near the critical probability obeys the following 
scaling relation (for p < Pc): 

5(p,/) = 5(pe,0-exp(^). (21) 

Where Apc = y. Given a system with a maximal chemical length /, for all values of 
conductivity p inside the range [pc — ApcPc + Apc] the survivability behaves similar to the 
power law S{pc, I) ~ found at criticality. Thus, the width of the critical threshold is 
Ap, = s^. 

To summarize, we have shown analytically and numerically the the survivability in ER 
and SF graphs scales according to equations (|TT|l and (fr2|) near the critical point. This 
implies that the width of the critical region in networks of finite size scales as Apc = y, 
where / is the chemical length of the percolation cluster. For ER graphs, / ~ N^^^, while for 
SF networks with 3 < A < 4, / ~ n{x~3)/{\-i) _ 
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FIG. 1: A graphical sketch of the recursion relation 0. Given a certain conduction probability p, 
the probability to reach i nodes at layer / + 1 is found by summing up the probabilities to follow a 
link to a node (whose outgoing degree is described by Gi{x)) and reaching i nodes at layer I from 
that node. 

University Press, to appear, 2005). 
[19] R. Cohen, D. ben Avraham, and S. Havlin, Phys. Rev. E 66, 036113 (2002). 
[20] O. Riordan and P. L. Krapivsky, private communication. 

[21] Indeed, solving equation (|10|) taking 5 <C 1 and the initial condition €^=0 = 1 we get: ei = 



[22] This method assumes that near the critical point there is a negligible number of loops and 
thus the network behaves similar to a cayley tree with the same degree distribution P{k) as 
the ER network. 
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site is conducting, multiplied by the probability Poo that it resides in the giant component. 

[24] Poo obeys the transcendental equation: Pqq = 1 — 6"^^^^^^°°. 

[25] Solving equation (^0]) with 6 <^ 1 and the initial condition e;=o = 1 we get: ei = 





[26] In this range is the behavior is similar to ER gi 
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FIG. 2: (Color online) The survivability S{p, I) for an ER graphs with {k) = 5, numerically 
calculated for different values of p: Pc, Pc±5- 10-^ Pc ± 3 • 10-^ Pc±l- 10""^, Pc ± 6.66 • 10-^ and 
j)c±3.33- 10^^. For p = Pc the survivability decays to zero according to a power law: S{pc, I) ~ 
For p < Pc, S{p, I) — > 0, while for p > pc, S{p, I) Const. The decay is exponential (to zero or to 
a constant) according to equations (fTT|) and dJ. 
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FIG. 3: (Color online) Scaling of the survivability for different values of p, I, and {k). Shown is 



S{p,l)-S{p,oc) 



VS. —\p — Pc\l for ER graphs with {k) = 5 (unfilled symbols) and (k) = 10 (filled 



symbols). The collapse of all curves on an exponential function (for large /) shows that indeed the 
scaling relations and ()12() are correct. 



10 



10 40 160 

I 



FIG. 4: (Color online) The survivability S{p, I) for a SF network with A = 3.5, numerically 
calculated for different values of p: pc, Pc i 6 • 10~^, pc± 4: ■ 10~^, pc ±2 ■ 10~^, 1.33 ■ 10~^, and 
Pci6.66- 10~^. For p = Pc the survivability decays to zero according to a power law: S{pc, I) ~ 
For p < Pc, S{p, I) — > 0, while for p > pc, S{p, I) Const. The decay is exponential (to zero or to 
a constant) according to equations (fTT|) and dJ. 
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FIG. 5: (Color online) Scaling of the survivability for different values of p, I, and A. Shown 
S{p,i)^^S{p,oo) — Pc\l for SF graphs with A = 3.5 (filled symbols) and A = 5 (unfilled 

symbols). Due to numerical difficulties only curves with p < pc are shown All curves collapse 
on an exponential function according to relation 

''The exact enumeration method is limited here to small chemical distances I due to the upper cutoff of 
the scale-free distribution. 
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